Spin polarized current and Andreev transmission in planar 
superconducting-ferromagnetic Nb-Ni junctions 
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We have measured and modeled, using three different approaches, the tunneling current in a 
Nb/Nb^Oy/Ni planar tunnel junction. The experimental data could be fitted and the correct 
current polarization could be extracted using a simple quasiclassical model, even in the absence of 
an applied magnetic field. We also discuss the microscopic structure of the barrier. 

PACS numbers: 74.50.+r, 74.80.Fp, 75.70.-i 



I. INTRODUCTION 

Experiments with tunnel junctions using ferromagnetic 
metals— have been an interesting topic since a long time. 
This subject has grown again 2 , because of the new field 
of spintronic where spin dependent currents are an im- 
portant requisite of many possible devices^*^ This im- 
plies that control and measurements of spin polarized 
currents are needed. Spin-polarized electron tunneling 5 
is a key tool to measure the current polarization and to 
understand the physics involved in these effects. Most of 
the recent experimental works have been focused on the 
suppression of the Andreev reflection and point contact 
geometry using a superconducting electrodei&lHowever 
the local information extracted by point contact or scan- 
ning tunneling microscope techniques seems to be less 
suitable for devices than planar tunneling junctions. In 
addition, the intrinsic difficulty of fabricating a perfect 
uniform oxide layer can jeopardize the latter technique. 
Recently Kim and Moodera— have reported a large spin 
polarization of 0.25 from polycrystalline and epitaxial Ni 
(111) films using Meservey and Tedrows's technique 5 , and 
standard Al electrode and oxide barriers. 

In this work we show that the tunneling electron po- 
larization of ferromagnets could be extracted without ap- 
plying a magnetic field to the junction and using an ox- 
ide barrier with random metallic point contacts. The 
experimental data are obtained for Nb/Nb^Oy/Ni pla- 
nar tunnel junctions. We analyze possible models which 
describe the experimental results, with emphasis on the 
information that can be obtained about the shape of the 
barrier, and the relation between the polarization of the 
transmitted current and the bulk polarization. 



II. EXPERIMENTAL METHOD 

Nb(110) and Ni(lll) films, grown by dc magnetron 
sputtering, were used as electrodes. The structural char- 
acterization of these films was done by x-ray diffraction 
(XRD) and atomic force microscopy (AFM), see for in- 



stance Villegas et al? Briefly, the junction fabrication 
was as follows: First, a Nb thin film of 100 nm thick- 
ness was evaporated on a Si substrate at room temper- 
ature. An Ar pressure of 1 mTorr was kept during the 
deposition. Under these conditions, the roughness of the 
Nb film, extracted from XRD and AFM, is less than 0.3 
nm— and superconducting critical temperatures of 8.6 K 
are obtained. After this, the film was chemically etched 
to make a strip of 1 mm width. A tunnel barrier was 
prepared by oxidizing this Nb electrode in a saturated 
water vapour atmosphere at room temperature,- Fig. ^ 
shows the x-ray diffraction data (open circles) of a Nb ox- 
idized film. The thickness of the oxide layer was deduced 
after the simulation of the data (solid lines) using the 
SUPREX program—i It can be seen that the oxide film 
thickness is 2.5 nm. X-ray photoelectron spectroscopy 
analysis performed in these oxidized films reveals that di- 
electric Nb205 is the main oxide formed, as can be seen 
in Fig. |3 There are also other oxides, such as metal- 
lic NbO, but in much less amount. Taken into account 
Grundner and Halbritter studies^ Nb205 is the outer- 
most oxide layer on Nb, whereas NbO is located closer 
to the Nb film. The characterization by AFM reveals a 
RMS roughness of around 0.7 nm. 

On top of this film (Nb with the oxide barrier), the 
second electrode of Ni was deposited under the same con- 
ditions as Nb (up to 60 nm thickness) using a mask to 
produce cross strips of 0.5 mm width, so that the overlap 
area S of the two electrodes is 0.5 mm 2 . 

Perpendicular transport (tunneling configuration) was 
investigated by means of characteristic dynamic resis- 
tance (dV/dl) versus voltage (V) using a conventional 
bridge with the four-probe method and lock-in tech- 
niques. The measured lock-in output voltage was cali- 
brated in terms of resistance by using a known standard 
resistor. 

Experimental data (open circles) in Figs. 3, 4 and 
6 show the behavior of the normalized conductance 
G(V)/G n versus the bias voltage V. Conductance has 
been calculated as the inverse of the differential resistance 
dV/dl and the normalization has been done with respect 
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FIG. 1: Measured (open circles) and simulated (solid line) 
low-angle x-ray diffraction profiles of a Nb film oxidized in 
a saturated water vapour atmosphere. The thickness of the 
oxide layer is 2.5 nm, as extracted from the simulation. 




200 



214 212 210 208 206 204 202 200 198 
Binding Energy (eV) 



FIG. 2: X-ray photoelectron spectroscopy analysis of a Nb 
film oxidized in a saturated water vapour atmosphere. Notice 
that the highest contribution component oxide comes from 
Nb 2 O s . 



to the dependent voltage background conductance G n 



III. THEORETICAL DESCRIPTION 

A. Introduction 

Conductance G across a normal-superconducting junc- 
tion may be expressed as a function of applied voltage V 
in terms of the reflection coefficient R of electrons travers- 



ing the junction as 
G{V) = 2 e 2 vp N F S I de (1 + A - R) 



df(e - eV) 
de 

where /(e) is the Fermi function; vf and Np, the Fermi 
velocity and density of states at the Fermi level, respec- 
tively; and A is the coefficient for Andreev reflection pro- 
cesses whereby an electron with energy smaller than the 
gap impinging onto the junction picks an electron of op- 
posite spin to form a Cooper pair inside the supercon- 
ductor, leaving behind a hole. Andreev reflection pro- 
cesses are proportional to the square of the conventional 
transmission coefficient of the barrier T, and therefore 
are strongly supressed for highly resistive barriers. Junc- 
tions with transmission coefficients smaller than about 
0.1 show small subgap conductances. 

Our experimental results, shown in Fig. [21 exhibit a 
significant conductance below the superconducting gap 
even at the lowest temperatures. Therefore we expect 
that our effective oxide barriers should be neither too 
high nor too thick. 

We have used three models to describe the transmis- 
sion across the oxide barrier: i) The simple generaliza- 
tion of the Blonder-Tinkham-Klapwijk (BTK) modeli^ 
to ferromagnetic electrodes proposed by Strijkers and 
coworkers^ ii) a generalization of the BTK model to 
ferromagnetic electrodes with finite bulk magnetization; 
and iii) a description of the effects of a finite current 
polarization in terms of spin dependent transmission co- 
efficients, in a similar way as discussed by Perez- Willard 
et alM.. 



B. Strijkers' model (Model I) 

Strijkers' model uses as adjustable parameters the cur- 
rent polarization 



Pr 



G T (0) + G X (0)' 



(2) 



the height of the barrier Z, which is modeled by a delta 
function, and the size of the superconducting gap at the 
interface A. The process of electron transfer across the 
junction is split into a fully polarized channel, for which 
the Andreev reflection is zero, and a paramagnetic chan- 
nel described by the BTK model in its usual form. The 
total conductance is then written in terms of the conduc- 
tance of the fully polarized channel (Gp) and the con- 
ductance of the paramagnetic channel (Gn)' 



G (V) = (1 - P c ) G N (V) + P c G P (n 



(3) 



This model interpolates between the paramagnetic case 
(BTK model), and the half metal, where it predicts cor- 
rectly that the amplitude for Andreev reflection vanishes. 

The best fit obtained with this model is shown in Fig. 
The parameters used are A = 1.4 meV (supercon- 
ducting gap), Z = 1.15, and P c = 0.03. This value of the 
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FIG. 3: Normalized tunneling conductance as a function of 
voltage of a Nb-Ni planar junction. Symbols correspond to 
experimental data and continuos solid line to the calculated 
curve using Model I. The temperature of the experiment was 
1.52 K and the fitting parameters were Z = 1.15, P c = 0.03, 
A = 1.4 meV (see text for details). 



polarization, however, is well below the usually reported 
current polarization of Ni. 



C. Generalization of BTK model for a 
ferromagnetic electrode (Model II) 

We next introduce ferromagnctism through 
an exchange splitting J in one of the elec- 
trodes. Therefore wave-vectors depend on spin as 
hk a — (2m(Ep + a J/2)) 1 / 2 . Transmission across 
the barrier is therefore spin-dependent, and Andreev 
reflection mixes both spins, giving different normal and 
Andreev reflection probabilities for each spin. Conduc- 
tance can be calculated for each spin in the same way as 
in the BTK model using these reflection probabilities. 
Total conductance will be then G = Gf + Gi, where 



G a = 2 e 2 v Fa N Fa S 



de(l + ^A B 



Ra) 



df(e - eV) 
de 



(4) 

The adjustable parameters in this model are J, the height 
of the barrier for zero splitting Z, again modeled by a 
delta function, and the gap at the interface A. 

This model leads to the results shown in Fig. 0] A good 
fit is obtained using Z = 1.13, A = 1.4 meV and J = 0.8 
eV, which agrees well with the electronic band structure 
of Ni. This exchange provides a bulk polarization 



J 



n-\ + ni 2 Ep 



0.4 



(5) 



if we take Ep ~ 1 eV, which is the botton of the band of 
Ni along our experimental TL direction. The model also 



FIG. 4: Fit to the experimental results obtained with model 
II (see text for details). The fitting parameters were Z — 1.13, 
A = 1.4 meV, J = 0.8 eV and the transmission coefficients 
were T T = 0.48 and T x = 0.37. 



allows us to define the current polarization, assuming a 
unidimensional barrier, in terms of the difference in the 
transmission coefficients through the barrier, Tf = 0.48 
and T i = 0.37. We find 



Pr 



0.13, 



(6) 



which we still feel to be somewhat too low. 

We note that the polarization as determined from the 
transmission coefficients needs not agree with the bulk 
polarization of the ferromagnetic electrodeiiS^i Similar 
effects can arise when the conduction bands are built up 
of localized and delocalized atomic orbitalsii^ A compar- 
ison between the polarization defined using Eq.jnj and 
the bulk polarization, for barrier strenghts in the range 
from 0.2 to 3.5 is shown in Fig. |SJ In the limits Pb — > 
and Pb — > 1 both polarizations are equal, but for interme- 
diate values current polarization is always smaller than 
bulk one. 



D. Simple quasiclassical theory (Model III) 

We now use quasiclassical theory, with boundary con- 
ditions dumped into the transmission coefficients T a . 
Conductance is then determined using Eq. J3J with 
fcf = ki, where the effective reflection coefficients are 
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1 + r a r_ 
-a + - r- a ) Q 



(1 - TV r_ 
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f + r CT r- a + (1 
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FIG. 5: Polarization determined from the transmission co- 
efficients across the barrier, eq.@ as a function of the bulk 
polarization of the ferromagnetic electrode. 



FIG. 6: Fit to the experimental results obtained with model 
III (see text for details). The fitting parameters were A = 
1.38 meV and the used transmission coefficients were Tj =0.5 
and T y = 0.32. 



Here g and / are the (spin independent) normal and 
anomalous components of the Green's functions evalu- 
ated right at the interface, at the superconducting side, 
and r 2 a +T a = 1. 

This scheme is, in principle, exact, although it does not 
allow us to extract information about the nature of the 
barrier, which is treated as a black box whose internal 
structure is ignored. We also assume that the electrodes 
are in the clean limit, to simplify the calculation. Once 
the transmission coefficients across the barrier which best 
fit the experimental data have been determined, we have 
performed calculations in the normal regime with more 
realistic square barriers in order to check the physical 
properties of the model. Neglecting effects related to the 
dispersion in the direction parallel to the junction, the 
results obtained using this model should be reproduced 
by the generalized BTK model with the appropiate choice 
of transmission coefficients. 

The best results obtained using this model are shown 
in Fig. H3 The parameters used are A = 1.38 meV, 
Tf = 0.5 and Tj_ = 0.32. Using the expression in eq.® 
we obtain a current polarization P c = 0.22. 



of about Ep « 1.5 eV along the TP direction. We find 
transmission coefficients comparable to those required to 
fit the experimental data for barrier depths d ~ 5 A. 

The significant conductance observed at voltages be- 
low the gap implies the existence of good contacts. One 
may get an estimate of the number of these contacts by 
dividing the quantum unit of resistance, h/(2e 2 ) by the 
observed resistance, R = 20 fi (the fact that the trans- 
mission coefficients needed to fit the data are below 1 
does not alter the order of magnitude). Assuming that 
the area of each of these conducting channels is of the 
order of a few tens of A 2 we obtain that the area of these 
contacts is a small fraction of the area of the oxide barrier 
(10 -8 ), suggesting that the conductance is mainly due to 
tiny spots where the barrier is narrower than the average 
thickness of the oxide layer. The existence of zones where 
the barrier is ~ 5 A thick is not inconsistent with the ex- 
perimentally observed corrugation, and agrees with our 
estimation of the barrier depth. 



IV. SUMMARY AND CONCLUSIONS 



E. Characterization of the barrier 

We have also related the transmission coefficients used 
in the fits with the thickness d and height U of more 
realistic square barriers. Grundner and Halbritter have 
found that the effective height of the barrier of NbaOs 
oxide layers is of order U w 0.1 eVi& Fig. [7| shows the 
transmission through a barrier of such height as function 
of d, for different values of bulk polarization. The Fermi 
energy in the Ni ferromagnetic electrode is approximately 
Ep w 1 eV, while in the Nb superconducting electrode is 



The experimental results can be fit with different mod- 
els. The simplest one, which describes the junction in 
terms of a paramagnetic and a fully polarized channel, re- 
quires as an input a bulk polarization which is much lower 
than the observed polarization of Ni. A more accurate 
one is the exact generalization of the BTK model, which 
describes properly the bulk polarization of Ni but not its 
current one. A scattering approach, which does not re- 
quire knowledge of the structure of the barrier, describes 
well the experimental data with reasonable values for the 
transmission coefficients for majority and minority spins, 
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FIG. 7: Transmission coefficients as a function of barrier 
depth for a barrier of heigth U — 0.1 eV. Bulk polarization 
in the ferromagnetic electrode is, from left to right and from 
top to bottom: 0.2, 0.4, 0.6, 0.8, as indicated in each graph. 



giving the correct value of the current polarization. The 



experimental data show an additional contribution which 
is not symmetric with respect to the applied voltage, this 
small asymmetry is not relevant for our results. This con- 
tribution may arise from inelastic scattering of magnons 
at the interfaced 

The polarization inferred from the transmission coeffi- 
cients calculated with the BTK model does not coincide 
with the bulk magnetization of the ferromagnetic elec- 
trode, used as an input. We have studied the relation 
between these two quantities, showing a significant de- 
pendence on the height of the barrier. 
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